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ABSTRACT 

The radiation released at the transparency radius of an ultrarelativistic flow can account for the 
observed properties of gamma-ray bursts (GRBs) provided that sufficient energy is dissipated 
in the sub-photospheric region. Here, I investigate how the peak energy of the Ef(E) spectrum 
and its overall shape depend on the properties of the jet for various "dissipative photospheres". 
I find that continuous energy release which results in electron heating over a wide range 
of distances may be the key to explain the GRB emission. In this picture, the peak of the 
spectrum forms at a Thomson optical depth of several tens. The peak depends mainly on 
the bulk Lorentz factor T of the flow and can, therefore, be used to determine it. The T is 
predicted to range from ~ 10 to 1000 from X-ray flashes to the brightest observed GRBs in 
agreement with recent observational inferences. The Amati relation can be understood if the 
brightest bursts are the least baryon loaded ones. Implications from this interpretation of the 
GRB emission for the central engine are discussed. 
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1 INTRODUCTION 

Although several thousand gamma-ray bursts (GRBs) have been 
observed so far, the mechanisms responsible for the GRB emission 
remain elusive. Both synchrotron and photospheric models have 
been widely explored in the literature to explain the characteris- 
tic "Band-like" GRB spectrum which is characterized by smoothly 
connected power-laws (Band et al. 1993). The synchrotron models 
rely on a power-law electron distribution accelerated at large dis- 
tance (or small optical depths) in the jet. The photospheric models 
focus, instead, on radiation that is released when the jet becomes 
transparent. 

Photospheric models are attractive interpretation because they 
can result in high radiative efficiency and naturally predict peak 
energies £ pca k ~ 1 MeV close to the observed values (Goodman 
1986; Thompson 1994; Meszaros & Rees 2000). The value of £ pea k 
in this interpretation is coupled to the main properties of the flow 
(e.g. luminosity L, and Lorentz factor T). Indeed, observations in- 
dicate that the peak of the Ef{E) spectrum tracks the instantaneous 
gamma-ray luminosity and integrated energy during a burst and 
among different bursts, respectively (e.g., Amati et al. 2002; Yo- 
netoku et al. 2004; Ghirlanda, Nava & Ghisellini 2010). 

On the other hand, the energy dissipated in the base of the jet 
effectively thermalizes, so in the absence of additional dissipation 
at modest optical depths (i.e. further out) the emission from the 
transparency radius is quasi-thermal in sharp contrast to that typi- 
cally observed. Significant dissipation of energy of some sort is re- 
quired close to the photosphere of the flow to lead to the observed, 
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smoothly-connected power-law spectra. The source of such dissi- 
pation may be (strong or many, weak) shocks (Rees & Meszaros 
1994; Lazzati & Begelman 2011), magnetic reconnection (Gian- 
nios 2006), or nuclear collisions (Beloborodov 2010). 

Adopting the reconnection model for GRBs of Drenkhahn 
(2002), I have shown that magnetic dissipation leads to powerful 
photospheric emission (Giannios 2006). The observed ~MeV peak 
of the spectrum forms at Thomson optical depth of several tens 
where radiation and electrons drop out of thermal equilibrium; the 
electrons turn hotter further out in the flow. Inverse Compton scat- 
tering at larger distance (smaller optical depth) leads to the high- 
energy tail that can extend well into the GeV range. A Band-like 
spectrum naturally forms in this scenario. 

Here, I develop a generic dissipative photospheric model ap- 
plicable to arbitrary dissipative process that results in electron heat- 
ing. Both a flat rate of dissipation of energy (e.g. constant lumi- 
nosity dissipated per decade of distance) and localized dissipation 
events are explored (Sections 2, 3). The model predicts a relation of 
the peak £ pea k of the observed emission with the properties of the 
flow; most sensitively depending on the bulk Lorentz factor T (Sec- 



tion 4). Using the observed E, 



peak 



L v relation I make inferences 



for the central engine of GRBs (Section 4.1). Section 5 clarifies 
various aspects of the GRB variability in the context of the model. 
Discussion and conclusions are presented in Section 6. 



2 DISSIPATIVE PHOTOSPHERES 

Energy dissipated at the very inner parts of the jet flow quickly ther- 
malizes. A substantial thermal component can also be built in an 
initially rather cold, magnetically dominated flow when magnetic 
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energy is dissipated at large Thomson optical depth t » 1, e.g. 
as expected in a striped wind model (Drenkhahn & Spruit 2002). 
The thermal luminosity L t \, (dominated by radiation) is released at 
the transparency radius (defined as the distance at which r = 1). 
For typical parameters of the jet flow, the resulting quasi-thermal 
emission peaks at E pcak ~ 3k^T obs ~ 0.1-1 MeV (Goodman 
1986, Thompson 1994; Daigne & Mochkovitch 2002; Beloborodov 
2010), very close to the typically observed values (Band 1993). In 
the absence of dissipation of energy close to the photosphere, how- 
ever, the emerging emission cannot account for the observed GRB 
spectrum. Though isolated cases for a strong quasi-thermal com- 
ponent in the GRB emission have been made (Ryde 2005; Ryde 
et al. 2010; Ryde et al. 2011), the GRB spectrum genetically has 
non-thermal appearance. 

The photospheric emission is, however, modified when addi- 
tional energy release takes place close to the transparency radius. 
It turns out (see next Section) that continuous electron heating at 
a range of optical depths from t ~several tens out to ~0.1 may be 
the key to reproduce the observed emission. Continuous dissipation 
results in a well defined distance where radiation and particles drop 
out of equilibrium, the so-called equilibrium radius. The peak of 
the emission spectrum is determined by the plasma temperature at 
this distance. Below we develop a general framework to calculate 
the peak of the spectrum as function of the properties of the flow. 



2.1 A generic model 

Consider a jet coasting with bulk Lorentz factor T, total isotropic 
equivalent luminosity L and baryon loading 77 = L/Mc 2 > T. In the 
presence of energetic particles injected by the dissipative process, 
the flow can be loaded with a modest number of pairs (Pe'er et 
al. 2006; Vurm et al. 2011). Assuming f ± electron+positron pairs 
per proton (i.e. /+ = 1 for e - p plasma; hereafter pairs are re- 
ferred to as electron^, the rest-frame electron number density is 
n e = f ± L/47tr 2 nTm p c 3 . Clearly both T and /+ can vary with distance 
(because of acceleration of the flow and energetic particle injection 
that result in pair creation, respectively). Here, I treat these quanti- 
ties as constants with the main focus been on their values close to 
the equilibrium radius. The Thomson optical depth t = n e cr T r/T as 
function of distance is 
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where all quantities are in cgs units and the A = A* 10* notation is 
adopted. Setting r = 1 , the Thomson photosphere is located at 
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We consider a flow which carries a thermal component of lu- 
minosity L t h that is a substantial fraction e of that of the flow: L t h = 



(4/3)4nr 2 r 2 aTl 



eL. This can be realized in both fireballs and 



dissipative, magnetically dominated flows (e.g. Drenkhahn 2002) 
The (rest frame) plasma temperature is 



knT, h 
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Inspection of eqs. (2) and (3) reveals that the plasma typically cools 
to sub-keV temperature close to the photosphere. It turns out, how- 
ever, that the dissipative process heats up the flow to a temperature 
well in excess of before it reaches the photosphere. 

Suppose that a dissipative process injects energy in the flow 
heating the electrons (see Section 3 for discussion on the physical 
justification of such assumption). I consider two different situations 
for the radial dependence of the dissipation rate: i) a gradual rate 
of energy release of rather flat profile (i.e. constant rate of energy 
dissipated per e-folding of distance): dL A jdr = L Ao jr and ii) dissi- 
pation that is localized to a narrow range in distance. 

In the case of the gradual energy release, the heating rate per 
unit volume in the rest frame of the flow is = L^al^nTr 1 . If the 
thermal component L th is built by the dissipation of energy from 
smaller distances (and including adiabatic cooling of the photons) 
one gets L,h = (3/2)Ld_ . Allowing for a comparable amount of 
additional heating even deeper in the flow, we set L Ao = = 
eL/3. 

Heating of the electrons is balanced by radiative cooling. 
Electrons and photons are in thermal equilibrium at large depth 
with their common temperature given by eq. (30. This equilibrium 
is broken once heating and cooling balance leads to T e > T th . It 
can be shown (see Giannios 2006) that Compton cooling Pic = 
4« c c ccr x f/ r is the dominant cooling mechanism for the electrons 
in this region (0 C = k^T c lm e c 2 and U T is the energy density of ra- 
diation). Equating dissipative heating and cooling (P^ = P\c) and 
setting U T = aT* h one derives the electron temperature as function 
of distance 
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The location where radiation and electrons drop out of equilibrium 
is found by setting = T s , defining the (maximum) equilibrium 



radius r„ 
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At the equilibrium radius the Thomson optical depth and 
the temperature of the electrons are, respectively 
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Note that the optical depth and temperature of the plasma at 
the equilibrium radius depend very weakly on the flow parameters 
and are in the range of several tens and ~ 1 keV, respectively. The 
equilibrium radius is of particular importance since it is the location 
at which the peak of the spectrum forms under a wide range of 
conditions. The observed peak of the spectrum is 



j peak 
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MeV. 
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2 Adiabatic cooling can be shown to be negligible for the electrons 
In reality, the electrons maintai 

see Ioka et al. (201 1). ation because of external heating 



1 For a photospheric GRB model where the flow does not contain baryons 3 In reality, the electrons maintain a slightly higher temperature from radi- 
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At distance r > r cq , and using eqs. (1) and (4), the Compton y 
parameters found to be y - 40 c t ~ 0.4 independently of distance 
or the parameters of the flow leading to significant up-scattering 
of the photons. The Compton parameter is close to unity because 
the incoming radiative luminosity and dissipation rate Ld, are 
comparable. Up-scattering of the radiation emerging from r cq at 
larger distance leads to broader spectra and to a flat high-energy 
tail above £ pea k. The high-energy tail is followed by an exponential 
cutoff at the energy that corresponds to the electron temperature at 
the distance where the dissipation stops (see next section). 

2.2 Numerical Results 

The analytical arguments presented in the previous Section are fully 
supported by detailed Monte-Carlo radiative transfer simulations. 
The code developed in Giannios (2006) is used to study the electro- 
magnetic spectrum emerging from a dissipative photosphere. The 
inner boundary of the calculation is set at the equilibrium radius 
where photons following a thermal distribution of temperature T cq 
are injected. The photons are followed throughout the photosphere 
until the flow reaches an optical depth of r = 0.1 where the outer 
boundary is set. The calculation includes Compton scattering, rela- 
tivistic effects, while the electron temperature is iterated until the 
heating rate is matched by Compton cooling everywhere in the 
flow. We do not include synchrotron emission and non-thermal par- 
ticle acceleration. While both are important in determining the ex- 
act spectrum, they are model dependent. For this paper, however, 
it is important that the peak of the emission is set at the distance 
where radiation and matter drop out of equilibrium and rather inde- 
pendently of the details of non-thermal processes. In the following, 
both extended dissipation with dLj/dr = h^ jr and localized dis- 
sipation are investigated. 

2.2.1 Extended Dissipation 

Assuming extended dissipation of energy, in Fig. 1 I plot the ra- 
diation spectrum at various distances of the flow for the reference 
values of the parameters (L = 10 53 erg-s -1 , r] = T = 300, e = 0.3, 
ft = 1). Note that the thermal emission at r eq evolves into one of 
non-thermal appearance when passing through the t = 1 surface 
building a flat high-energy tail. Inverse Compton scattering with 
Compton y parameter y < 1 results in a flat (Ef(E) ~ E°) emis- 
sion above the peak. The outer boundary of the calculation is set at 
larger distance that corresponds to r = 0. 1 . The high-energy cut- 
off E cat is determined by the temperature of the electrons at the 
outer boundary (in this example kgT^t ^ 300 keV resulting in 
£cut - r£ B r collt ^ 100 MeV). While in this paper, I set the outer 
boundary of dissipation by hand at r = 0.1, more detailed models, 
such as the magnetic reconnection model of Drenkhahn & Spruit 
(2002) predict where this cutoff takes place. 

In Fig. 2 I show the resulting spectra for different values of the 
parameters. The peak of the emission E pcak is very close to the ana- 
lytic expression (8). The overall shape of the spectra is very similar 
with only the peak and cutoff energies depending on the parame- 
ters. A Band-like spectrum is reproduced for very different param- 
eters of the flow. The high-energy tail has spectral slope f(E) oc EP 
with p < - 1 as observed. Below the peak energy the slope, as mea- 
sured in the 10-100 keV range, is f(E) oc E a with a ~ 0-1. 

4 For relativistically expanding fluid the number of scatterings scales as ~ r 
and not oc t 2 as in a static medium (see Giannios 2006). 
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Figure 1. Numerically calculated spectra at different optical depths in the 
jet for the reference values of the parameters and for extended dissipation 
of energy. The thermal injected photon spectrum at the equilibrium radius 
(black line) evolves into a broader spectrum at the r = 1 surface (red line) 
and develops a flat power-law tail at the outer radius (corresponding to t = 
0.1; green line). The peak of the Ef(E) emerging spectrum is determined at 
the equilibrium radius. 




E (keV) 



Figure 2. Emerging spectra for different values of the parameters and for 
extended dissipation of energy. The spectrum has a Band-like shape practi- 
cally independently of the luminosity, baryon loading and Lorentz factor of 
the flow. The peak of the spectrum follows the scaling predicted by eq. (8). 
The high energy cutoff is set by the location of the outer boundary of the 
simulation. 

The slope a is consistent but in the rather hard range in compari- 
son with the observed distribution. Note, however, that synchrotron 
and synchrotron-self-Compton emission from larger distances (not 
included in this work) can soften the spectrum below the peak (Gi- 
annios 2008; Vurm, Beloborodov & Poutanen 201 1). 

2.2.2 Localized Dissipation 

The process that dissipates energy in the jet may result in energy 
release over a narrow range of distances. Here, I investigate how a 
thermal component carried by the flow is modified depending on 
the optical depth at which such energy release takes place. 

I assume that the thermal component L,h = eh is reprocessed 
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E (keV) 

Figure 3. Emerging spectrum for dissipation of energy at a narrow range 
of distance (or optical depth T&ss) for different Tdi 8s . Dissipation at large 
optical depths r t |j ss > 10 leads to narrow emission spectrum. Dissipation at 
Tdiss S; 1 leads to distinct thermal and high-energy components. For r&a, ~ 3 
the two components are smoothly connected. 

by Compton up-scattering through a narrow region of hot electrons. 
The electrons are heated at a rate L i o = L th over a range of Thom- 
son depth T diss ...T diss /2 (i.e. a factor 2 in distance). The electron tem- 
perature is determined by heating-cooling balance. Setting the var- 
ious parameters to their reference values, Fig. 3 shows the emerg- 
ing spectrum for different values of Tdj ss . Dissipation at large opti- 
cal depths Tdj ss > 10 leads to a broadened, Planck-like distribution. 
The resulting narrow emission spectrum has a steep rise/decline 
below/above the peak. Dissipation at T^ is < 1 leads to a hot re- 
gion above the photosphere. Most of the photons do not interact 
with the hot electrons leading to a distinct quasi-thermal compo- 
nent while; photons that are up-scattered at least once form a high- 
energy component. Although occasionally GRB spectra show such 
multi-component behavior in the keV-MeV regime, this in not typ- 
ical. For Tdi K ~ 3 the two components are smoothly connected and 
the resulting spectrum compares more favorable with observations. 
For a model where the dissipation (through magnetic reconnection) 
takes place at a narrow region of t > 1 see Thompson (1994). 

Summarizing, localized dissipation in general does not ac- 
count for observations because it either leads to narrow or multi- 
component spectra in the ~MeV energy range. There is a limited 
range of optical depths of Tdi SS ~ 3 - 5 where dissipation results in 
a smooth Band-like spectrum. 



3 THE DISSIPATIVE MECHANISMS 

In the previous analysis we simplified the calculation assuming 
thermal electrons continuously heated by an external agent. Here 
we elaborate on the possible sources for smooth "volume" heat- 
ing of the GRB flow (Ghisellini & Celotti 1999; Stern & Poutanen 
2004; Pe'er, Meszaros & Rees 2006). 

Models for gradual energy release that heats the electrons 
involve magnetic reconnection (Giannios 2006), multiple weak 
shocks (Ioka et al. 2007; Lazzati & Begelman 20100 and neutron- 

5 The, more often invoked for the GRB emission, strong internal shocks 
(Rees & Meszaros 1994) are unlikely to lead to smooth/gradual heating of 
the electrons. Such shocks, instead, accelerate particles practically instan- 



proton collisions (Beloborodov 2010). In the first case, the energy 
is initially stored in the magnetic field while in the latter cases in 
relative bulk motions within the jet. In the following I argue that, 
independently of the dissipative mechanism, a large fraction of the 
energy dissipated close to the equilibrium radius is expected to be 
stored to (mildly) relativistic protons. Coulomb e — p collisions ef- 
fectively drive the energy from protons into the electrons (see Be- 
loborodov 2010). The electrons then transfer their energy into the 
photon field through Compton scattering. The characteristic elec- 
tron equilibrium temperature is sub-relativistic and determined by 
heating-cooling balance. 

Multiple weak shocks in the jet or elastic neutron-proton col- 
lisions (close to the distance where the neutron and proton fluids 
decouple) can be expected to heat the protons at mildly relativistic 
temperatures with the heating maintained over a range of distances 
in the jet. When the heating takes place at optical depth of tens, 
the Coulomb coupling of < GeV protons with ~keV electrons (see, 
e.g., Stepney 1983) can be shown to be effective in transferring the 
energy into the electrons within one expansion timescale of the jet. 
These, slowly heated, electrons pass their energy into the photon 
field effectively through inverse Compton scattering. In addition to 
mildly relativistic protons, shocks and inelastic n — p collisions can 
inject high-energy particles. Non-thermal processes can result in 
a modest pair loading in the jet (Pe'er et al. 2006; Beloborodov 
201 of] but not expected to directly affect the peak energy of the 
emission as long they are not energetically dominant. 

Magnetic reconnection can take place effectively at large op- 
tical depths (deep in the jet) through, e.g., tearing instabilities of 
the current sheet (Loureiro, Schekochihin & Cowley 2007; Uzden- 
sky, Loureiro & Schekochihin 2010) assisted/induced by the accel- 
eration of the jet (Lyubarsky 201 1)0 Magnetic reconnection takes 
place in multiple dissipative centers that accelerate particles, heat 
the plasma and drive fast bulk motions. The fast motions in the 
downstream of the reconnection region are required for sufficiently 
fast reconnection. The bulk motions are dissipated further down- 
stream in shocks resulting in hot protons. Coulomb collisions can 
effectively couple the proton energy to the electrons which, in turn, 
couple it to the radiation field (as discussed above). In this strongly 
magnetized plasma (t/ B > (7 ph ), sub-relativistic electrons can ther- 
malize very efficiently through exchange of synchrotron photons 
(the so-called "synchrotron boiler"; Ghisellini et al. 1998) before 
they have the chance to cool through inverse Compton scattering. 
An alternative picture where the energy released in the reconnec- 
tion process drives MHD waves (instead of bulk motions) is de- 
scribed in Thompson (1994). Photons extract the energy from the 
waves through scattering on electrons (Compton drag). Also in this 
picture the electron velocity is limited to sub-relativistic speed and 
close to a Compton equilibrium with the radiation field. The analy- 
sis of the previous section may apply to this scenario as well. 

In summary the details of the dissipation process (and the re- 
lated non-thermal processes) may well affect the pair injection rate 
and result in additional features in the spectrum (e.g. powering a 
pair annihilation line at E ~ Tm t c 2 ) but do not change the basic 

taneously to ultrarelativistic energies leaving them to cool radiatively on a 
longer time scale (see, however, Ramirez-Ruiz 2005). 

6 Note, however, that if the energy density of the magnetic field is higher 
than that of the radiation, ultrarelativistic electrons cool mainly through syn- 
chrotron emission and the pair creation is further limited. 

7 A possible speed up of the reconnection rate at Thomson thin conditions 
(Lyutikov & Blandford 2003; McKinney & Uzdensky 2012) can have inter- 
esting implications but is not required in this picture. 



© 201 1 RAS, MNRAS 000,[TJ-?? 



Peak energy of dissipative GRB photospheres 5 



picture of sub-relativistic electrons strongly coupled to the radia- 
tion field with their energy determined by a balance of dissipative 
heating and radiative cooling. It is this basic process at the equilib- 
rium radius that determines the peak of the emerging emission. 



4 INFERRING THE BULK LORENTZ FACTOR OF THE 
FLOW 

In this rather general framework of dissipative photospheres, there 
is a close connection between the peak of the spectrum and the 
properties of the flow summarized in eq. (8). Eq. (8) can be solved 
in terms of T: 

r-osns- 3 ' 5 -1/10, 1/10 f 1/5/ 1/5 /q\ 
1 - 280£ pcakMcV e_ 05 L 53 f ± {-J , (9) 

where I chose to express T as function of the ratio n/T > 1 instead 
of r\. A clear implication from the last expression is that T depends 
mainly on the peak energy E w ^ with extremely weak dependence 
on the physical properties of the jet. To first approximation, one 
can have a fair estimate of the Lorentz factor just from the observed 
peak of the GRB emission. From eq. (9) it is clear that the model 
predicts that the high-peaked GRBs come from faster jets. This is 
in sharp contrast to the synchrotron internal shock model where the 
opposite holds true (£ pca k K T~ 2 ). 

One can somewhat improve on the estimate of the bulk T of the 
jet from observables by assuming a gamma-ray luminosity L y ~ eL 
and using the observed L y and e ~ 1. This interpretation of the 
peak of the emission implies that the brightest GRBs with Epeak ~ 
several MeV and L 53 ~ 10 (e.g. Abdo et al. 2009) come from the 
most relativistic T ~ 1000 flows while weaker X-ray flashes with 
£peak ~ 30 keV and L53 ~ 0.01 come from "slower" jets of T ~ 20. 
Low-luminosity GRBs (e.g. Soderberg et al. 2006) and the X-ray 
flares that follow many bursts (Burrows et al. 2005; Chincarini et 
al. 2010) can also be a result of a dissipative photosphere from yet 
slower jets of T < 10. In particular low-luminosity GRBs may come 
form T ~ 3 jets that result in only modest relativistic beaming of 
their emission and may, therefore, account for the larger observed 
local rate of these events (Soderberg et al. 2006) in comparison to 
classical GRBs. As discussed in Section 6, observational estimates 
of the Lorentz factor of GRBs support the model prediction that 
higher £ pea k bursts are coming from higher T flows. 



4.1 Additional inferences for the flow using observed 
correlations of the GRB emission 

It has been recognized for some time that various time integrated 
quantities over the duration of a burst may correlate with each- 
other, e.g., the peak energy £ pca k with the isotropic gamma-ray en- 
ergy E y (Amati et al. 2002). Even in the presence of outliers (Band 
& Preece 2005; Nakar & Piran 2005), these correlations may teach 
us a lot about the GRB physics. The emission in the photospheric 
models is, however, connected to the instantaneous properties of 
the flow. A change in any of the properties (e.g. luminosity or 
baryon loading) during the burst shifts the location of the photo- 
sphere and its appearance. Ghirlanda et al. (2010; 201 1) found that 
there is a time dependent Amati-like relation of £ pc ak(f) and L y (t) 
during the evolution of bursts where E^.^ ^ MeV. 

If this relation is verified by more data, it implies for the pho- 
tospheric models that the Lorentz factor correlates with the lumi- 
nosity of the flow. In this interpretation, one can derive an estimate 



of the Lorentz factor of the flow directly from the observed gamma- 
ray luminosity using eq. (9), the instantaneous E^^-Ly correlation 
and that L y ~ L lh = eL: 

r = 310<i, McV e-^'^)" /5 . (10) 

Note that the Lorentz factor depends sensitively on a single observ- 
able, namely, the peak energy. Therefore, E vcilk can be use the infer 
the Lorentz factor. In terms of properties of the flow, I find that 
T = 200e'_Q 5 L 2 5 ^(r]/ry l/5 fl 15 , i.e., that more luminous bursts are 
less baryon loaded with T oc L 2/s . Note, however, that a systematic 
dependence of any other quantity with, say, the luminosity of the 
flow can distort this T - L scaling. Finally, the T - L relation can 
be combined with eq. (7) to find the temperature of the flow at the 
equilibrium radius: 

k B T eq = O.S^L^ /;"^)" 5 keV. (11) 

Note that the temperature of the flow clusters in the ~keV range 
and has very weak dependence on the parameters of the jet. As 
discussed in Section 6 such clustering of the peak of the emission 
in the rest frame of the flow has been observed. 

The instantaneous E ps . A - L y relation also puts interesting con- 
straints on the distance where the peak of the spectrum forms. Us- 
ing the expression (5) for r cq and eq. l |10t , I arrive to 

r eq = 2.1xl0"L-^ 1 Vr/i /5 (^[ 1/5 cm. (12) 

The equilibrium distance r cq depends very weakly on the various 
parameters. Even allowing for orders of magnitude variations in 
the luminosity and (as expected) more modest changes in other pa- 
rameters from burst to burst (or during the evolution of a burst), r eq 
varies at most by a factor of a few. 



5 GRB VARIABILITY 

GRBs are variable on timescales as short as milliseconds. Their 
lightcurves are characterized by multiple pulses (e.g. Fenimore 
et al. 1995) that typically last for seconds and show characteris- 
tic spectrum that often evolves from hard to soft (e.g. Hakkila & 
Preece 2011) and/or tracks the instantaneous flux of the flow (e.g. 
Ghirlanda et al. 201 1). Can the photospheric model account for the 
temporal GRB behaviour? 

The small radius of emission in the photospheric models 
allows for fast variability down to f v ~ r ph /2r 2 c ~sub-msec 
timescales (depending on the parameters of the flow; see also Gian- 
nios & Spruit 2007). Allowing the continuous heating out to optical 
depth t ~ 0. 1 (r out = 10r ph ) can lead to the multi MeV high-energy 
tail delayed by r out /2r 2 c ~several msec with respect to the MeV 
peak (with the GeV emission potentially more delayed depending 
on the distance at which it takes place). The steady-state jet model 
developed here is not applicable to study extremely short (sub- 
msec) timescales (for which the steady-state assumption breaks 
down) but is well suited to study the evolution of the burst on longer 
timescales. In particular, the ~sec duration GRB pulses can be ac- 
curately studied as a sequence of steady state models where the 
properties of the jet (L = L(t), F = T(t), etc) vary on this timescale. 
The observed variability on second timescales in the well studied 
~ 10- 1000 keV energy range is likely to be dominated by temporal 
evolution of the properties of the jet rather than delays introduced 
by propagation effects of the ejecta. 

The spectral evolution during a GRB pulse depends on how 
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T(t) and L(t) evolve during the pulse. Giannios & Spruit (2007) 
have shown that If, for instance, T(t) has a positive power-law de- 
pendence on L(t) during individual pulses then the peak of the spec- 
trum tracks the observed flux (as seen observationally in Ghirlanda 
et al. 201 1). In this case, while the flux declines after the peak lumi- 
nosity of a pulse, the peak of the emission spectrum also declines 
and the spectrum softens. As a result the pulse duration on the softer 
X-ray bands is longer than in the harder ones. 

In summary, continuous dissipation close to the photosphere 
can allow for fast evolving (sub-msec) emission. Furthermore, it is 
possible to account for observed temporal properties of GRBs given 
specific assumptions about the behaviour of the central engine. The 
reasons for which the engine operates in such a fashion is, however, 
not addressed by this work. 



6 DISCUSSION AND CONCLUSIONS 

In this paper, I have explored a wide range of "dissipative photo- 
sphere" models to pin down the location where the peak iipeak of 
the spectrum forms and how it connects to the properties of the 
jet. The resulting expression (9) of my analysis summarizes this 
connection, is quite general and rather independent of the physical 
model for energy dissipation. 

I have explored both localized and continuous in distance dis- 
sipation of energy. Localized heating of electrons over a narrow 
range of distance is shown here to have difficulties to account for 
observations. If the dissipation takes place at large optical depth, 
the emerging spectrum is rather narrow, while dissipation at r < 1 
has two distinct components in the ~MeV energy range, both in 
conflict to the majority of observed bursts. However, localized dis- 
sipation that takes place at optical depth t ~ 3 results is a more 
promising spectrum. Thompson (1994) discusses a scenario where 
such dissipation is possible. 

I find that continuous dissipation of energy over a wide range 
of optical depths can naturally give the Band-like spectrum with 
peak and slopes above and below the peak at the observed range. 
Magnetic reconnection in a jet that contains small-scale field rever- 
sals naturally results in such flat dissipation profile throughout the 
photospheric region. Alternatives such as multiple, weak shocks or 
neutron-proton collisions can also result in continuous energy in- 
jection. The MeV peak of the spectrum forms at the distance where 
radiation and the electrons drop out of thermal equilibrium. In the 
context of this model the peak energy is mainly determined by the 
bulk Lorentz factor of the flow. The observed E p!:ak can, therefore, 
be used to infer T. 

The model predicts that the peak energy £ pca k positively cor- 
relates with T. I find that weak GRBs (the spectrum of which peaks 
in the X-ray regime; the so called X-ray flashes) and X-ray flares 
that follow the bursts potentially come from T ~ 10 jets while the 
brightest observed Fermi-LAT bursts, with peak energy at several 
MeV, come from the fastest T ~ 1000 flows. Recently (Ghirlanda 
et al. 2012) used the peak time of the afterglow emission to estimate 
the bulk Lorentz factor T of 31 GRBs. They verified that brighter 
bursts are characterized by higher T and showed that the peak en- 
ergy of the emission at the rest frame of the jet clusters at several 
keV; both findings are unique predictions of this work (see eqs. 7, 
11). 

The observed relation of the peak of the spectrum with 
the instantaneous burst luminosity indicates a close coupling of 
the emerging spectrum to instantaneous properties of the flow 
(Ghirlanda et al. 2011). It has already been pointed out in Gian- 



nios & Spruit (2007) that the - L y relation implies, in the con- 
text of the reconnection model, that brighter segments of a burst 
come from higher T (i.e., "the brighter the cleaner") jet. I come 
to the same conclusion for the generic photospheric model devel- 
oped here where the bulk Lorentz factor T and the luminosity of 
the flow L scale as, roughly, T oc L 2/5 . Note again that independent 
constraints coming from afterglow modeling indicate a similar cor- 
relation for T and the gamma-ray energy E y (Liang et al. 2010). 

In this model, the MeV peak forms in a fairly compact region, 
while the high-energy tail forms at a larger distance. The tail can 
extend to multi-GeV energy without suffering attenuation due to 
pair creation and can exhibit delays of the order of r Ge v/r 2 c with 
respect to the MeV emission that may be of order of seconds (as 
observed in some Fermi-LAT bursts, see Abdo et al. 2009; ro e v 
is the radius where GeV emission takes place). While a different 
mechanism may be invoked for the dissipation at large distance that 
leads to the GeV emission (Meszaros & Rees 201 1), the continuous 
spectrum from MeV to GeV energy indicates a single dissipative 
mechanisrrf] (see, e.g., Bosnjak & Kumar 2012). 

Interestingly, using the observed instantaneous £ pea k - L y re- 
lation, the equilibrium distance r eq (where the peak of the emission 
is set) can be shown to depend extremely weakly on the parameters 
of the flow with r tq = 2 x 10" cm. This value is similar to the radii 
of Wolf-Rayet stars, probable progenitors of long-duration GRBs. 
This could indicate that interactions of the jet with the progeni- 
tor are important for energy dissipation (e.g. Lazzati & Begelman 
2010). On the other hand, a similar correlation is observed in short- 
duration GRBs, supporting the idea of a progenitor-independent 
dissipative mechanism (Ghirlanda et al. 201 1). 

The inferred "the brighter the cleaner" property of GRB flows 
has profound implications for the nature of the central engine. Both 
accreting black holes and rapidly rotating magnetars have been in- 
voked for launching the jet. Models that invoke accretion into a 
black hole (Woosley 1993) are not sufficiently developed to predict 
the amount of baryons that make it into the jet. In the millisecond 
magnetar model for GRBs (Usov 1992; Metzger et al. 2011), the 
calculation of the baryon loading is more tractable. The magnetars 
born with the strongest fields drive the brightest bursts and also 
give, averaged over the GRB duration, less baryon loaded flows. In 
particular, oblique rotators result in T oc L° 6 (but with a large scat- 
ter; see Metzger et al. 2011). This is a rather intriguing result since 
it might point to a rather complete picture for GRB physics (which 
is, however, hardly unique). A central engine consists of a proto- 
magnetar (with its magnetic field axis, in general, misaligned to the 
rotational axis), that gives rise to a magnetically dominated jet that 
contains field reversals on small scale (striped wind). Magnetic re- 
connection proceeds in the jet at a wide range of distances from the 
central engine and at a rather flat rate (Drenkhahn & Spruit 2002). 
The bright MeV emission of the GRB emerges from the Thom- 
son photosphere of the flow while residual dissipation can extend 
the high-energy tail above the MeV peak well into the GeV regime 
(Giannios 2008). 
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8 but not necessarily a single emitting region for the MeV and GeV emis- 
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